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Abstract 

We consider sensor transmission power control for state estimation, using a Bayesian inference approach. A sensor node sends 
its local state estimate to a remote estimator over an unreliable wireless communication channel with random data packet 
drops. As related to packet dropout rate, transmission power is chosen by the sensor based on the relative importance of 
the local state estimate. The proposed power controller is proved to preserve Gaussianity of local estimate innovation, which 
enables us to obtain a closed-form solution of the expected state estimation error covariance. Comparisons with alternative 
non-data-driven controllers demonstrate performance improvement using our approach. 
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1 Introduction 

Wireless networked systems have a wide spectrum of ap¬ 
plications in smart grid, environment monitoring, intel¬ 
ligent transportation, etc. State estimation is a key en¬ 
abling technology where the sensor(s) and the estimator 
communicate over a wireless network. Energy conserva¬ 
tion is a crucial issue as most wireless sensors use on¬ 
board batteries which are difficult to replace and typ¬ 
ically are expected to work for years without replace¬ 
ment. Thus power control becomes crucial. In this work, 
we consider sensor transmission power control for re¬ 
mote state estimation over a packet-dropping network. 
Transmission power control in state estimation scenario 
has been considered from different perspectives. Some 
works took transmission costs as constant. Shi et al. [1] 
assumed sensors to have two energy modes, allowing it 
to send data to a remote estimator over an unreliable 
channel either using a high or low transmission power 
level. The optimal power controller is to minimize the 
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expected terminal estimation error at the remote esti¬ 
mator subject to an energy constraint. Similar works can 
also be found in [2,3]. Meanwhile, some literature has 
taken channel conditions into account. Quevedo et al. [4] 
studied state estimation over fading channels. They pro¬ 
posed a predictive control algorithm, where power and 
cookbooks are determined in an online fashion based 
on the undergoing estimation error covariance and the 
channel gain predictions. More related works can been 
seen in [5-7]. 

An important issue which has not been taken seriously in 
most works is that the transmission power assignment, 
as a tool to control the accessibility of information to the 
receiver, should be determined not only by the underly¬ 
ing channel condition and the desired estimation perfor¬ 
mance, but also by the transmitted information itself. 
In [4] and [5], the authors failed to associate transmis¬ 
sion power with data to be sent. The plant states are 
used to determine the transmission power in [8]. In this 
case, lost packets signal the receiver of the state infor¬ 
mation. To avoid computation difficulty, the signaling 
information is discarded. 

In this paper, we focus on how to adapt the transmis¬ 
sion power to the measurements of plant state and how 
to exploit information contained in the lost packets. We 
propose a data-driven power controller, which utilizes 


Preprint submitted to Automatica 


25 March 2015 



different transmission power levels to send the local esti¬ 
mate according to a quadratic function of a key param¬ 
eter called “incremental innovation” which is evaluated 
by the sensor at each time slot. By doing this, even when 
data dropouts occur, the remote estimator can utilize 
the additional signaling information to refine the pos¬ 
terior probability density of the estimation error by a 
Bayesian inference technique (see [9]), therefore deriving 
the MMSE estimate. It compensates the deteriorated es¬ 
timation performance caused by packet losses. To facili¬ 
tate analysis, we assume that a baseline power controller 
has already been established based on different factors 
with regard to different settings, such as the requirement 
of estimation performance as in [1] or the channel condi¬ 
tions as in [4,5,7]. We are devoted to developing a power 
controller that embellishes this baseline controller by 
adapting the transmission power to the measurements 
such that the averaged power with respect to all possible 
values taken by the measurements does not exceed that 
of the baseline power controller. The proposed power 
controller, driven by online measurements, can run on 
top of non-data-driven power controllers, which results 
in hierarchical power control mechanisms. Then exten¬ 
sion to a time-varying power baseline is established in 
Section 4.4. Note that a related controller was first pro¬ 
posed in [10], but as a special case of the controller in 
this work. The main contributions of the present work 
are summarized as follows. 

(1) We propose a data-driven power control strat¬ 
egy for state estimation with packet losses, which 
adapts the transmission power to the measured 
plant states. 

(2) We prove that the proposed power controller pre¬ 
serves Gaussianity of the local innovation. It sim¬ 
plifies derivation of the MMSE estimate and leads 
to a closed-form expression of the expected state 
estimation error covariance. 

(3) We present a tuning method for parameter design. 
Despite of its sub-optimality, the controller is shown 
to perform not worse than an alternative non-data- 
driven one. 

The remainder of this paper is organized as follows. 
In Sections 2 and 3, we give mathematical models of 
the considered system and introduce the data-driven 
transmission power controller. In Section 4, we present 
the MMSE estimate at the remote estimator and a 
sub-optimal power controller that minimizes an upper 
bound of the remote estimation error. In Section 5, com¬ 
parisons with alternative non-data-driven controllers 
demonstrate performance improvement using our ap¬ 
proach. Section 6 presents concluding remarks. 

Notation: N (and N+) is the set of nonnegative (and 
positive) integers. is the cone of n by n positive 
semi-definite matrices. Eor a matrix X, A^(X) is the ith 
smallest nonzero eigenvalue. We abuse notations det(X) 
and which are used, in case of a singular matrix 

X, to denote the pseudo-determinant and the Moore- 


Penrose pseudoinverse. Sij is the Dirac delta function, 
i.e., Sij equals to 1 when i = j and 0 otherwise. The no¬ 
tation pdf(x, x) represents the probability density func¬ 
tion (pdf) of a random variable x taking value at x. 

2 State Estimation using a Smart Sensor 

Consider a linear time-invariant (LTI) system: 

Xk+i=Axk^Wk, ( 1 ) 

Vk ~ Cxk T Vk-, (2) 

where /c G N, G is the system state vector at time 
k, yk ^ is the measurement obtained by the sensor, 
the state noise Wk G and observation noise Vk G 
are zero-mean i.i.d. Gaussian noises with E[wkWj] = 

SkjQ {Q h 0), E[i;/e('z;j)'] = SkjR {R >- 0), E[wk{vjy] = 

0 Vjf, /c G N. The initial state xq is a zero-mean Gaus¬ 
sian random vector with covariance Bq ^ 0 and is un¬ 
correlated with Wk and Vk. (B, C) is assumed to be de¬ 
tectable and (B, Q^/^) is assumed to be stabilizable. Eur- 
t her more, we assume A is Hurtwitzp~] 



Fig. 1. The system architecture. 

2.1 Sensor Local Estimate 

Hovareshti et al [11] illustrated that utilization of the 
computation capabilities of wireless sensors may im¬ 
prove the system performance significantly. Equipped 
with such “smart sensors”, the sensor locally runs a 
Kalman filter to produce the MMSE estimate x^. of the 
state Xk based on all the measurements collected up 
to time k, i.e., yi^k = {^i,..., ^/c}, and then transmits 
its local estimate to the remote estimator. Denote the 
sensor’s local MMSE state estimate, the corresponding 


^ Since we focus on remote state estimation in this paper, 
for any practically working systems (to be monitored alone), 
A has to be Hurwitz. Otherwise, the system state will go 
unbounded and there is no real sensing device which can 
track an unbounded state trajectory. Adding a control input 
to regulate the system state for an unstable A and studying 
its associated stability issue will be beyond the scope of this 
paper and will be left as our future work. 
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estimation error and error covariance as and 

respectively, i.e., xl = E[xk\yi:k], el = Xk - xl and 
= E[{xk — xl){xk — xiy\yi:k]‘ Standard Kalman 
filtering analysis suggests that these quantities can be 
calculated recursively (cf., [ 12 ]), where the recursion 
starts from xg = 0 and Pq = Ho ^ 0. Since P^ converges 
to a steady-state value exponentially fast (cf., [ 12 ]), we 
assume that the sensor’s local Kalman filter has entered 
the steady state, that is, P^ = P P 0 V/c G N, This as¬ 
sumption simplifies our subsequent analysis and results, 
such as Theorem 4.8 and Proposition 4.17. 

2.2 Wireless Communieation Model 

The data are sent to the remote estimator over an Addi¬ 
tive White Gaussian Noise (AWGN) channel using the 
Quadrature Amplitude Modulation (QAM) whereby xl 
is quantized into K bits and mapped to one of 2^ avail¬ 
able QAM symbols p~| For simplicity, the following as¬ 
sumptions are made: 

A.l: The channel noise is independent of Wk and Vk. 

A. 2 : P is large enough so that quantization effect is neg¬ 
ligible when analyzing the performance of the re¬ 
mote estimator. 

A.3: The remote estimator can detect symbol errorj"^ 
Only the data arriving error-free are regarded as 
being successfully received; otherwise they are re¬ 
garded as dropout. 

These assumptions are commonly used in communica¬ 
tion and control theories (cf., [4, 5, 8 , 13,14]). For ex¬ 
ample, Fu and Souza [14] demonstrated that the esti¬ 
mation quality improvement (in terms of reduction of 
the remote estimation error) achieved by increasing the 
number K of the quantization bits is marginal when K 
is sufficiently large (in their example K only needs to 
be greater or equal to 4. Based on A.3, the communica¬ 
tion channel can be characterized by a random process 
{7fe}feeN+, where 

_ Jl, if arrives error-free at time /c, 

\o, otherwise, 

initialized with 70 = 1. Denote ji^k — { 71 , • • • ,7/c}- Let 
uJk G [0, -hoc) be the transmission power for the QAM 
symbol at time k. We adopt the wireless communication 
channel model used in [10], and have Pr ( 7 /^ = 0|cc;/c) = 
^ where q is given hy q = exp(—(a/(Poll^)) ^ ( 0 ? 1 )^ 
Nq is the AWGN noise power spectral density, W is the 
channel bandwidth, and a G (0,1] is a constant that de¬ 
pends on the specific modulation being used. To send lo¬ 
cal estimates to the remote estimator, the sensor chooses 


^ QAM is a common modulation scheme widely used in 
IEEE 802.11g/n as well as 3G and LTE systems, due to its 
high bandwidth efficiency. 

^ In practice, symbol errors can be detected via a cyclic 
redundancy check (CRC) code. 


from a continuum of available power levels cj/c ^ 0 , see 
Fig. 1 . Note that different power levels lead to different 
dropout rates, thereby affecting estimation performance. 

2.3 Remote State Estimation 

Define Ik as the information available to the remote es¬ 
timator up to time k, i.e., 

h = (3) 

Denote Xk and Pk as the remote estimator’s own MMSE 
state estimate and the corresponding estimation error 
covariance, i.e.^Xk =E[xk\lk] andP/c = E[{xk—Xk){xk — 
XkYl^k]^ where expectations are taken with respect to a 
fixed power controller. We assume that the remote esti¬ 
mator feedbacks acknowledgements yk before time kpl. 
Such setups are common especially when the remote es¬ 
timator (gateway) is an energy-abundant device. This 
energy asymmetry allows the estimator to trade energy 
cost for estimation accuracy. 

3 Data-driven Transmission Power Control 

Our strategy uses the measurements to assign transmis¬ 
sion power level efficiently. As focusing on how the power 
controller utilize the sensor’s real-time data, to simplify 
discussion, we assume a constant power baseline 00 in 
this section. We define 0 = {0k}keN+ as a transmission 
power controller over the entire time horizon, where Ok 
is a mapping from yi^k and yi^k to uJk. Before proceeding 
to study 6 >, let us first briefly explain the idea of data- 
driven power control mechanism. Define r(/c) G N+ as 
the holding time since the most recent time when the 
remote estimator received the data from the sensor, i.e., 

T{k) = k- : 7t = 1}. (4) 

We interchange r(/c) with r when the underlying time 
index is clear from the context. Define Sk as the incre¬ 
mental innovation in the sensor local state estimate com¬ 
pared to time k — T^ the previous reception instant, i.e., 

£k=xl- (5) 


Lemma 3.1 E[e|e^|Ifc_i, 7 fe = 0] = 0 Vfc G N+. 

Proof: The result follows from noting that 

E[efeefe|Ife_i,7fc = 0] = E[E[ele'^.\yi:k,'ri-.k]\h-i,'rk = 0] 

= E[E[efe|2/i:fe]efc|Ifc-i,7fe = 0] = 0, 

where the seeond equality holds beeause el is independent 
of yi:ky ond the last equality holds sinee E[el\yi:k] = 0 . 
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Note that, if Sk = 0, then the sensor generates a lo¬ 
cal estimate, identical to the prediction We 

would say that, for the remote estimator, the “value” 
of information contained in x^. is null. As Sk becomes 
larger, x^ has an increasing drift from the prediction 
and the importance of the sensor sending x^, 
thereby raises. Motivated by these observations, we de¬ 
fine a stationary power controller, 6 >ef : 5 /c ^ cj/c, as 
an increasing function of Sk- To fit the above observa¬ 
tions, we introduce a quadratic function of Sk given by 
C(5/c, Q) = Sk Qsk^ where Q G is a weight matrix. 
According to Lemma 3.1, the covariance of Sk is a func¬ 
tion of r{k). Therefore we specify r(/c) for the index of 
Q and construct the following controller: 


or resorts to approximations, e.g., treating a truncated 
Gaussian distribution as a Gaussian distribution as was 
done in [15]. These approaches either lead to conserva¬ 
tive results (due to the unutilized information) or inac¬ 
curate results (due to approximations). Our method, on 
other hand, makes use of the information contained in 
the event 7 /^ = 0 to improve the estimation performance. 
The associated MMSE estimate, relying on no approxi¬ 
mation techniques, is derived in a closed-form. 

4 Main Results 

4.1 Preliminaries 


NnW 

6»ef; {ujk = —z — C{sk, Qt) + oj). (6) 

Za 

In contrast to ( 6 ), most non-data-driven transmission 
power controllers (i.e., [4,5]) use a given power uj regard¬ 
less of what value Sk takes. Note that in ( 6 ) a constant 
term uj is added after C{sk^ Qr)- If one sets = 0 , 
then the transmission with the baseline power controller 
cj = d) is a special case of the proposed transmission 
power controller. As for ^ 0 , the transmission power 
is a constant uj if C(5/cQr) = 0 ; otherwise it is adapted 
according to C{sk^ Qr)- Gompared with a related con¬ 
troller proposed earlier in [ 10 ], Oef in ( 6 ) is more gen¬ 
eral at least from two aspects: 1 ) we introduce a weight 
matrix Q^- to highlight the roles of different entries of 
Sk] 2 ) it allows the sensor to transmit using a stan¬ 
dard power UJ even if C{ek, Qr) = 0 , which includes a 
non-data-driven power transmission as a special case. 
As shown later in Lemma 4.4, given I/c_i, Sk is zero- 
mean Gaussian with a covariance depending on r(/c), 
he., {ek\fk-i) r\j A/’(0, S^-). For convenience of our subse¬ 
quent analysis, we define a new parameter satisfying 
Tt- = (^Qk + F“^) ^, where ^ ^ 0. We now list 

the main problems considered in the remainder of this 
work, 

(1) Under 6 >ef defined in ( 6 ), what is the MMSE es¬ 
timate and its associated estimation error covari¬ 
ance? 

( 2 ) What value should Qr (or T^-) take in order to 
minimize, 'E[Pk], the expected estimation error at 
the remote estimator? 

The solution to the first problem is presented in Sec¬ 
tion 4.2. A sub-optimal solution to the second one is 
given in Section 4.3 in view of the difficulty of the opti¬ 
mization problem. 


For any F ^ 0 that is singular, there exist matrices 
U^D G such that F = UDU'^ where U is uni¬ 

tary, whose columns are right eigenvectors of F, and 


D = 


"a 0 
. 0 0 


, where A is a diagonal matrix generated by 


the corresponding nonzero eigenvalues of F. Let F^/^ = 
U^/D. Then F = FV2 (fV2)' . 


Generally speaking, an n-dimensioned random vector 
X ^ A?(//, F), does not have a pdf with respect to the 
Lebesgue measure on if some entries in x degener¬ 
ate to almost surely constant random variables. To work 
with such vectors, one can instead consider Lebesgue 
measure in the rank(F)-dimension affine subspace: Q = 
{/r + F^/^z : z G with respect to which x has a 

pdf pdf(x, x) = ^ exp {—^{x — /i)'F“^(x — //)), where 

a = (27r)™^^^Met(F). Without loss of generality, in 
the remainder of this paper, for a random variable x ^ 
A/'(0, F) with a singular F, the pdf of x means the prob¬ 
ability density on Q. Note that the Moore-Penrose pseu¬ 
doinverse of F is unique and given by 



( 7 ) 


and that the pseudo-determinant of F equals to the prod¬ 
uct of all nonzero eigenvalues of F. 

Consider the power control law 6 >ef defined in ( 6 ). In or¬ 
der to guarantee that ujk is always nonnegative for any 
value 5/c, the difference of and F“^ needs to be 
at least positive semi-definite, i.e., two conditions must 
be simultaneously satisfied, which are F^- ^ and 
A F“^. The following lemma provides a necessary 
condition that needs to satisfy. 


Before proceeding, we note that in previous works such 
as [ 8 ] the difficulty of using the information contained 
in lost packets, i.e., 7 /^ = 0, when computing the MMSE 
estimate of the plant state has been acknowledged. One 
typically discards such information as was done in [ 8 ] 


Lemma 4.1 Supposed andT satisfyH P T A 

F“^. Then 

rank(T) = rank(F) ( 8 ) 

and 

Im(Ei/2) ^ (9) 
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where Im(X) is the image of X. 

Proof: Sinee , it is true ^/ia^rank(S) > rank(^). 

To verify {S), suppose that rank(i;) > rank(^). Then 
from (1), rank(S“^) > rank(^“^)^ whieh eontradiets 
with^!~^ ^ 11“^. To prove {9), let us denote idjak(^!) = r 
and assume there is a set of veetors W = {wi,..., w^} 
sueh that = span ({wi,..., w^}). Suppose 

Im(S^/^) 7 ^ Then there exists a veetor in W 

(without loss of generality, let it he mvi), and a veetor 
Wo G Ker where the operator Ker(X) is the 

kernel of a matrix X, sueh that wq'wi ^ 0. It leads to 
the faet that wq ^ Ker ((^^/^)'). IKe in turn have 

Wo = 0 while wo > 0, 


with ^0 = 0- dt 'Is also true that, given 7 /c = 0 and Ih-i, 
{ek\lk-u7k=0)^X{0,^r). 

Proposition 4.5 Under defined in (6), given 
Ik-i, the paeket drop rate at time k is given by 

Pr(7^ = 0|I*_,) = ^ 

We denote the packet arrival rate as — 1 ~ = 

0|I/c_i), where the subscript r is to emphasize that it 
depends on and 4^^-. To ensure that the averaged 
transmission power with respect to different values 
taken by the measurement in 9 does not exceed Co, i.e., 
E[co’/c|I/c-i] < Co, we require the following result. 

Lemma 4.6 Under (6), given Ik-i, the relation be¬ 
tween E[!jOk\lk-i] ci'o.d and uo is given by 


whieh eontradiets with S ^ T. 


E[a;fc|Ifc_i] = ^ (Tr(S,^'-i) - n,) +a;. (11) 


For convenience, denote — rank(ST-) = rank(TT-), 
^Iin(s7) =Im(^'7)and$^ = 

One has next lemma, the proof provided in the Ap¬ 
pendix. 

Lemma 4.2 The rank of^^ equals that (or 

i.e., rank(4>T-) = 

Example 4.3 Two matriees are provided below as a 
simple example for n = 3, 



1 

or 

O 

O 

_1 


’ 3 

-1 o' 

s, = 

0 5 0 

, and = 

-1 

3 0 


0 0 0 


0 

0 0 


Proof: From Lemma we know that {ek\lk-i) ^ 
A/’(0, S^-). Under Oq^, we have: 

E[LOk\h-i] = E [E[cj/c|5/c]|I/c_i] 


NqW 

2 a 

NqW 

2 a 

NoW 

2 a 


E[e', 

Tr(E[e,4|Ife_i] (v|>y - )) + 

(Tr(ET-5'“^) - rir) + w. 


00 


With 6>ef defined in (6), the remote estimator computes 
Xk and Pk according to the following two theorems. 


We can verify that Ur = 2, T,r h ^'r ^ ^ ^ (^)> Theorem 4.7 Under 0^i (6), the remote estimator com- 

and Lemma f.l holds. putes Xk as 


4.2 MMSE State Estimate 

In general, the posterior distribution of Sk fails to 
maintain Gaussianity without analog-amplitude obser¬ 
vations. The defect is especially common for quantized 
Kalman filtering and Gaussian filters, where it is tack¬ 
led by Gaussian approximation [12,16,17]. By contrast, 
the following lemma shows that, using ^ef in (b), the 
distribution of Sk conditioned on lk-i,^k = 0 is Gaus¬ 
sian. The proof, similar to that of Lemma 3.5 in [10], is 
omitted. 

Lemma 4.4 Under defined in (6), given Ik-i, ^k 
follows a Gaussian distribution: {ek\lk-i) ^ A/’(0, E^-), 
where is given by the following reeursion: 


xl, if 7fc = 1, 

A'^xl_.^, if 7fc = 0, 


( 12 ) 


where x^. is updated as x^. = -h Sk when 7 /^ = 1. 


Proof: When 7 /c = 1, the result is straightforward sinee 
xl is the MMSE estimate of Xk given yi:k. Now eonsider 
7 /c = 0. The tower rule gives 

E [a::fe|Ife-i,7fe = 0] = E [E [xhlvi-.k^i-.k] |Ifc-i,7fe = 0] 

= E [A'^xl_.^+ £fe|Ife-i,7fe = O] 

= ^^f7T + E[efe|Ife-i,7fe = 0]. 


S, = A^r-iA' + {h{P) - P) , 


( 10 ) 


Lemma 4-4 leads to E [5/c|I/c-i, 7/c = 0] = 0. 
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Theorem 4.8 Under (6), Pk at the remote estimator 
is updated as 


P, if 7k = 1, 

if 7k = 0. 


(13) 


Problem 4.10 


mm 




exp 


a 


-UJ 


NoW 


s.t. 


NqW 

2 a 


(Tr(ST-^'^ ^) - rir) + UJ < 


UJ. 


Proof: When 7 /c = 1 the result is straightforward. We 
only prove the ease when 7 /c = 0. 

^ [{a^k ^/c) {a^k ^/c) |l-/c— 1 5 7k ~ 0] 

= E [{Xk - A^xl_V{Xk - A^Xk_r)'\^k-l,lk = 0] 

= E [E [(e| + ek){-)'\yi-.k, li-.k] |Ifc-ij 7fc = 0] 

= E[(e|)(e|)'|yi.fc] +E[(efe)(£fe)'|Ife_i,7fc = 0] 

= P+^r, 

where the third equality is due to Lemma 3.1 and the last 
one is from Lemma 4-4- 

Remark 4.9 Under a baseline power eontroller with a 
eonstant power eontrol Co, the remote estimator’s esti¬ 
mate still obeys the reeursion (12); however, the estima¬ 
tion error eovarianee is updated differently: Pk = P when 
7 /e = 1, and Pk = h{Pk-i) = when 7 /e = 0^. Note 
that although the obtained estimates under the two power 
eontrollers are the same, their different estimation er¬ 
ror eovarianee matriees suggest different eonfident levels 
with whieh the remote estimator trusts the obtained es¬ 
timate: with the data-driven power eontroller, it is more 
eonvineed that the obtained estimate is elose to the real 
state while less eonvineed with a non-data-driven power 
eontroller. 


The constraint is imposed by (11). To solve Prob¬ 
lem 4.10, we first note that Tr(ST-T“^) = Tr(4>T-). 
However, for any matrix X,Y G det(XT) = 

det(X)det(T) does not hold in general since det(X) 
means X’s pseudo-determinant (in case X is singular). 
Fortunately, this property still holds for and 
The proof is given in the Appendix. 

Lemma 4.11 Suppose and satisfy Sr ^ ^ 0 

and y Then det(Sr)det(T:;r^) = det(<Fr). 

From linear algebra, det(<FT-) = 11^=1 

Tr(<F^) = We simply write Xi{^r) as 

Xi{r), and denote the nonzero eigenvalues of 4 >t- by 
Ar — [ Ai(r),..., (r) ]. Then Problem 4.10 can be 

recast as 

Problem 4.12 


min 

S.t. 


Ai(r)n”:iA.(r)V2 

NqW 


exp 


a 


-UJ 


2 a 


XI A* W 


NqW 
-\r UJ = UJ, C(J ^ 0 


(14) 


Li=l 

1 < Ai(r) < XjJ-), Vj = 2, 


4-3 Seleetion of Design Parameters 

The performances of ^ef foi* different T^-’s are difficult to 
compare in general. However, for Sr and Tr, there must 
exist a real number e^ G (0,1] such that Tr X ^rSr and 
Tt- eS^-, V e < er- Observe that 

= (s^y ^r-lEl/2 ^ 1 


lur 0 
0 0 


which yields er = 71(^7 • light of (10), we further 
have X ^rS^- = er + Si). According to 

Proposition 4.5, it can be seen given rfk) = r that 
E[P/c|t(/c) = r] has an upper bound: E[P/c|t(/c) = t] X 
P + (1 — p^)er (ATt-_iA' + Si) . Instead of minimizing 
E[P/c], we minimize its upper bound which is equivalent 
to minimize (1 — p^)er. Iterating over time, one eventu¬ 
ally needs to minimize (1 — p^)er for any T(k) G N+ at 
any k G N+. To this end, we propose to assign param¬ 
eters of 6>ef in (6) as the solution to the following opti¬ 
mization problem: 


Lemma 4.13 Let A* be the optimal solution to Prob¬ 
lem 4 U 2. Then Af satisfies 

Ai(r)*=A 2 (r)* = --- = A„^(r)*. (15) 


Proof: Suppose that A is the optimal solution to Prob¬ 
lem 4U2 but does not satisfy (15). IFe will show that 
there must exist another veetor, whieh is different from 
A and has a smaller eost funetion (14). Let Ai = c 
where c is a positive eonstant. Due to the faet that Ai 
in A is the minimum eigenvalue of 4>r and the inequal¬ 
ity of arithmetie and geometrie means, we have Ai < ^ 
/ \ 

and nr=i ^ equalities simultaneously 

satisfied when Xi = \f i = l,...,nr. Thus, Aq = 
..., results in a smaller value of (14), whieh eon- 
tradiets with the assumption and eompletes the proof. 


The following lemma is a result of Lemma 4.13. Its proof 
is presented in the Appendix. 
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Lemma 4.14 If UJ > , then the optimal solution to 

Problem 4-12 is uj = uj — and 

^ a 


A* — [1 H-,..., 1 H-]. 

Ur Ur 


(16) 


Otherwise, if uj < optimizer is uj = 0 and 


. ^ . 2auj ^ 2auj 

^r = [i + :r^7777---l + 


UrNoW 


UrNoW^ 


(17) 


determines the transmission power level, bit rates and 
codebooks used by the sensors. The algorithm in [4] re¬ 
quires that the receiver (i.e, the remote estimator) runs 
a channel gain predictor, see e.g., [18]. A key observa¬ 
tion is that the data-driven controller proposed in the 
present work can be readily adapted to situations where 
the baseline controller provides time-varying power lev¬ 
els u}/c|3]ln fact, by solving Problem 4.12 for a time- 
varying power baseline cD/c, we obtain the optimal solu¬ 
tion as follows: If ook > ? then is given by 


Denote by the transmission power associated with 
the solution to Problem 4.12. Then we have the following 
theorem. It can be readily verified from Lemma 4.14. 

Theorem 4.15 IfO> then 6>*£ is given by 


n* . _NoW , 

^ef • i^k — - ^k 


aUr 


_ NoW. 

-UJ -1, 

a 


where Hr+i = + h{P) — P with Hq = 0- 

Otherwise, if uj < is given by 

etf : {uj, = 

^here S.+i = + h(P) - P. 

Remark 4.16 A non-data-driven baseline power eon- 
troller with a eonstant power level uj is feasible to Prob¬ 
lem 400 . Sinee 0*^ is the optimal solution, it has not 
worse state estimation performanee eompared with the 
alternative non-data-driven power eontroller. Numeri- 
eal examples in Seetion 5 demonstrate performanee im¬ 
provements using 0*^ eompared with the non-data-driven 
power eontroller. 


^ef • 


{wfe = 




arij 




■ Wfe 


NqW, 


a 


and Otherwise, if uj, < 6'*f is 

given by 

Oli ■ K = 'eU 

and both cases, T., = (1 - 

7 /c_i)AT/c_iA' + h{P) — P. Note that H/., ^k and 
are calculated similar to and 4 >t- given in The¬ 

orem 4.15. To reduce the sensor’s computational load, 
the sensor only needs to calculate the quadratic form 
while the rest of the paraments are updated 
and then sent to the sensor by the estimator. Note that 
calculating has a complexity of O(n^). 


5 Simulation and Examples 


Consider a system with parameters as follows: A = 


0.99 0 . 3 ' 

,c = 

’ 2.3 1 

0.1 0.7_ 


. 1 


R = Q 


12 x2‘ We first 


assume that 0 has a constant power baseline u) = 5 and 
= 3 < In Section 5.2, a time-varying power 
baseline is considered. 


The following proposition shows that the rank of can 
be calculated offline. The proof is given in the Appendix. 

Proposition 4.17 Consider the given in Theo¬ 
rem 405, for any r G N+; ean be ealeulated as: 
Ur = rank(h'^(P) — P). In partieular, when r > n, the 
dimension of x, n^- beeomes a eonstant whieh is given 
by: Ur = rank(h’^(P) — P), \/ r > n. 

4.4 Extension 

In many cases, the base-line power controller changes 
over time with respect to different settings. For example, 
in [4], block fading channels were taken into account. To 
deal with a time-varying channel power gain h/Q a pre¬ 
dictive power control algorithm was established, which 


^ The term “channel power gain” means the square of the 
magnitude of the complex channel. 



Fig. 2. Empirical estimation covariance provided by con¬ 
trollers ^ef(^i) O 2 as a function of energy constraint Cd. 

5.1 Comparison with Different Energy Constraints 

We compare our proposed schedule 0*^ (denoted as ^ 1 ) 
with a constant baseline power controller within the en- 


^ Following assumptions commonly made in the literature, 
see, e.g., [4,7], in the sequel we shall assume that the chan¬ 
nel gain hk is available via the one-step ahead channel gain 
predictor. 
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30 


tire time horizon (denoted as O 2 : { 00 k = Define 

Jk{0) = ^ '^i=i Tr (IE[Pi]) as the empirical approxima¬ 
tion (via 100000 Monte Carlo simulations) of the aver¬ 
age expected state error covariance (denoted as 
We choose JsoiO) as an approximation of J{0). 

Fig. 2 shows that 6 >*£ leads to a better system perfor¬ 
mance when compared to O 2 under the same energy con¬ 
straint. 
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5.2 Comparison under Fading Channels 


Fig. 4. Comparison of ^ef(^i) O 3 given a specific realiza¬ 
tion of channel power gains. 


In practice, wireless communication channels typically 
comprise fading often assumed to be Rayleigh [19], i.e., 
the channel power gain hk is exponentially distribiRed 
with pdf(h/c) = iexp(—^), where hk ^ 0 and h is 
the mean of hk- Truncated channel inversion transmit 
power controllers have been studied in several works [5, 
7,20], where the transmission power is the inversion of 
h/c, with a truncated boundary. In this subsection, we 
use the baseline power determined by truncated channel 
gain inversion Denote the truncated channel inversion 
transmission power controller as ^ 3 : 


^k 


hk^ 

otherwise. 


(18) 


whnre v and h^ are design parameters. Consider the case 
of h = 1 and set h* = 5. Based on the results in [5], we 
can choose v to meet the energy constraint. Fig. 3 sug¬ 
gests that leads to better system performance when 
compared with ^ 3 . Fig. 4 shows the comparison given a 
specific realization of channel power gains. 



Fig. 3. Comparison of ^ef(^i) ^tnd O 3 under Rayleigh fading. 


6 Conclusion 

We proposed a data-driven transmission power con¬ 
troller for remote state estimation, which adjusts the 
sensor’s transmission power according to its real-time 
measurements. Then we proved that the proposed power 
controller preserves Gaussianity of the incremental in¬ 
novation and provided a closed-form expression of the 
expected state estimation error covariance, a tuning 


method for parameter design was presented to guar¬ 
antee that the data-driven power controller not worse 
performance than the alternative non-data-driven ones. 
Comparisons were conducted to illustrate estimation 
performance improvement. 


Appendix 


Proof of Lemma 4-2: To verify the clain, it suffices to 
show that rank(^T-) > Suppose that rank(^T-) = 
r < Ur. Since ^ 0 , there must exist exactly n — 
r mutually orthogonal vectors ei,...,en-r such that 
= 0, for i = 1,..., n — r. Denote the unit vec¬ 
tor with only the {nr5-j)th. entry being 1 by ij, that 
is, ij = [ 0 ,... , 0 , 1 , 0 ,... , 0 ]'. Since = 0 , j = 

’ v ^ 

rir+j 

1,..., n—without loss of generality, let = ij. As we 
assume that e^-r is orthogonal to ej, j = 1 ,..., n — n^-, 

it is true that D\^‘^en-r 7 ^ 0. Since Ur is nonsingu¬ 
lar and KeT{Ur) = { 0 }, we have e = 7 ^ 0 . 

We then observe that = en-r'^r^n-r = 0, and 


e'S, 


- 1 ^ 


^n — r 


In. 0 
0 0 


e^-r > 0 , which contradicts 


with ^ ^ H 


-1 


Proof of Lemma 4-11: By definition, it is easy to see that 
det(I]^)det(5'-i) = nr=iTherefore 

we only need to prove det($T-) = nr=i 

Observe that and can be factorized as = 


Ur 


0 

0 0 


Ur and Tt- 




0 ^ 0 

0 0 




where 


At- and 0 ^- are diagonal matrices generated respec¬ 
tively by the nonzero eigenvalues of and For 
i = 1 ,..., and Vi are the eigenvectors asso¬ 

ciated with A^(St-) and A^(Tt-). In addition, Ur = 


[ Ml,..., , 0 ,..., 0 ] and Vr = 




Then can be written as ^t- 


Mr 0 
0 0 


where Mr= 




























and Vr = ['L’l,...,'L’n^]- According to Lemma 4.3, 
is nonsingular, so det(<l>T-) = det(MT-). Since 
Im(ST-) = Im(4^T-) from (9), there exists a unitary 
matrix V such that Vr = UtV. Thus, det(MT-) = 
det j ^ (A^e^-^) , which 

completes the proof. ■ 


Proof of Lemma 4-14-' According to Lemma 4.13, we set 
Ai(r) = • • • = Xfirir) = At-. Logarithm does not change 
the monotonicity of (14). Problem 4.12 is consequently 
transformed to 


min 

Xt,UJ 

S.t. 


(X , Tir ^ \ 1 A 

c^-(- + l)lnA, 


rirNoW 

2a 


(At- — 1) (J — C(J, 


CJ > 0. 


(19) 


Substituting uj = + ^(At- — 1) — (^ + 1) In At- 

into (19) and taking derivative, it yields that the min¬ 
imum of (19) is attained at At- = 1 + Meanwhile 
uj needs to be nonnegative, so the optimal solution to 
Problem 4.10 is (16) if a) > or (17) otherwise. ■ 

Proof of Proposition 4jJ'^- Consider a matrix S = 
SI =1 Pi {h\P) - h*“^(P)) with pi G (0,1]. We have 


Im(I]) = ••• prA^-^T.\'‘^ ][■]’) 

= Im([E{/^ AE{/^ ••• A'^-'^t\^^]) 

= Im([Sy^ ••• ][■]') 

= Im(V(P)-P), (20) 


which leads to the first assertion. By the Cayley- 
Hamilton theorem, we have = —ai{k)A^~^ — 

a 2 {k)A^~‘^ - an{k)I, \f k > n, where ai(/c),..., an{k) 

are coefficients of the characteristic polynomial of A. 
When T > n + 1, we have 

Im([Ei/^ AE}/^ ••• ][•]') 

= Im([Ei/^ AE}/^ -ai(T-l)^”“iE}/^ 

-a2(T-l)A"-2E}/'-a„(r-l)Ep]), 

The last assertion follows from the reasoning used 
in (20). ■ 
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